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Hi ! 

We investigate the energy loss of quarks and gluons produced in hard processes resulting 
from final state rescatterings in a finite quark-gluon plasma. The angular distribution of 

7—i • the soft gluon bremsstrahlung induced by n s = 1 rescatterings in the plasma is computed 
in the Gyulassy-Wang model. Special focus is on how the interference between the initial 

t^j- ! hard radiation amplitude, the multiple induced Gunion-Bertsch radiation amplitudes, and 
gluon rescattering amplitudes modifies the classical parton cascade results. 
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Q\ • 1. INTRODUCTION 
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At RHIC energies (y/s ~ 200 AGeV) a new observable was predicted in nuclear col- 
lisions: jet quenching 0-fl. Collisional energy loss of a hard jet propagating through 
quark-gluon plasma has been estimated to be modest [|l],||: dE coU /dx <C 1 GeV/fm, 
however the induced radiative energy loss estimated HJU] from the Gunion-Bertsch am- 
plitude is expected to be significantly larger, dE rad /dx > few GeV/ fm. In Ref. a 
model for non-abelian energy loss was developed and the analog of the QED class of radi- 
ation diagrams was summed. This leads to a constant dE/dx ~ 1 GeV/fm independent of 
the path length. However, in BDMPS it was shown that gluon rescattering diagrams 
in the medium considerably increased the energy loss. For an incident very high energy 
jet penetrating a plasma of thickness, L, in which the mean free path, A = l/(ap), and 
the color electric fields are screened on a scale /z, the analytic BDMPS prediction || is 
dE/dx oc a s fi 2 L/\. An alternative path integral approach led to similar results 0. The 
transverse momentum dependence of the energy loss in QED was calculated recently |10| . 
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Unfortunately, the analytic approximations are not applicable to "thin" plasmas of only 
a few mean free path thickness, and more importantly they do not apply at all to the 
problem of computing the angular distribution. The main complication is that the angular 
distributions require a complete treatment of all gluon and jet final state interactions. In a 
direct brute force pQCD approach, in the thin plasma limit where the main simplicity 
arises from the fact that for a few rescatterings, n s < 3, the number of amplitudes is 
however still sufficiently small to allow direct computation of all interference terms. Here 
we display the complete investigation of n s < 3 rescattering can be found 

in Ref. |]]. 



The single rescattering case, n s = 1, was first considered by Gunion and Bertsch || 
(GB). They estimated the hadron rapidity distribution due to gluon bremsstrahlung as- 
sociated with a single elastic valence quark scattering in the Low-Nussinov model. The 
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most important difference compared to the problem addressed in this paper is that GB 
computed the non-abelian analog of the Bethe-Heitler formula while we are interested in 
the analog of radiative energy loss in sudden processes analogous to /3-decay. The GB 
problem therefore is to compute the soft radiation associated with a single scattering of 
an incident on-shell quark prepared in the remote past, i.e., t = — oo, relative to the 
collision time, ti. In our case, the radiation associated with a hard jet processes in nu- 
clear reactions must take into account the fact that the jet parton rescatters within a 
short time , t\ — to ~ Ra, after it is produced. The sudden appearance of the jet color 
dipole moment within a time interval ~ 1/E results in a broad angular distribution of 
gluons even with no final state scattering. The induced radiation cased by rescattering 
necessarily interferes strongly with this hard radiation. 

To model the scattering in a thin plasma, we employ the same model of the plasma 
as considered in Gyulassy-Wang (GW) H. The scattering centers are approximated by 
static color-screened potentials with Fourier components V(q) = 4:iia s / ((f + ji 2 ) and color 
screening mass, fi = A-Ka s T 2 . 

2. Angular Distribution with One Rescattering 

The amplitude for the production of a hard jet with momentum Pq localized initially 
near Xq = {to,0) we denote by M.j = J(P)e lPx °. The hard vertex is localized within 
a distance ~ 1/Po, and the amplitude, J(P), is assumed to vary slowly with P on the 
infrared screening scale /x. This jet radiates a gluon with momentum k, which can scatter 
on n s = 1 scattering center. 
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Fig. 1. Three contributions to the soft gluon radiation amplitude Mj® Mi. 

The conditional double inclusive jet and gluon probability distribution is the following: 

6 ( W .=l) _ >,=!),, , d*k d 3 p _ / |2\ ^ d*p 

d D rad - Prad (k, P ) (27r)3 ^ (27r)3 ^ -{\Mj® M ns=1 \ ) (27r)3 ^ (27r) 3 2 ^ ,(] 

which involves averaging over all initial target colors and summing over all final colors 
and the locations x« of the target centers. After averaging one obtains: 
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Tr (^^MUKp^^MljiKp;^ J . (2) 



The main source of non-diagonal dependence arises through the interference terms in 
and the T(qi±_ — q{ ± ) transverse form factor. 
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Following Ref. |§, the three radiative amplitudes Mu of Fig.l. are given by 



Mi tl = 2ig s ^^(e it ^-^)a 1 c, Mi, 2 = 2ig s ^-^{-e it ^)ca 1 , 

K i 



M h3 = 2tg s €± (u {K ^V 1 c x (e^^^ - e^^^) [c, oj . (3) 

The momentum distribution of gluon radiation can be derived from the conditional prob- 
ability distribution of radiation pf} d {p, k), normalized by the elastic scattering p^\p) fll!] . 
Before averaging on scattering centers one can obtain from eq. (|2|): 

i^ = P ^wrv- =► ^{^ 2 + ^(5 1 2 + t(5 1 )c ? 1 2 )-^(^-4cos(t 10 ^ )) 

-Rt(b x /2) (H ■ Ci cos(t w cu 10 ) - 2C l ■ Bi cosfaoWi)) j(4) 

Here we rearranged the radiation amplitudes and introduced the hard term if = k±/kj_, 
the cascading term Ci = (k — qi)±/(k — qi)\, and the GB term B\ = H — C\. The 
interference terms contain phase factors composed from transverse energy transfer terms 
ujq = k\/2u, lji = (k — qi) 2 j_/2u, u w = uj\ — oj§ and time difference term t w = t\ — t . 
These interference terms are responsible for "quantum cascading" : without them one can 
recover the classical cascading limit in eq. @; with them a complex behavior of "quantum 
cascading" dominated by these interference terms appears. The transverse form factor 
t(b) = T(6)/T(0) connects to the target size and it regulates automatically the gluon 
cascading on the finite size target object. 

In the limit t\ — > +oo the hard radiation formula for the no-scattering case, H 2 is 
recovered. The limit to ~~ * — oo is corresponding to the isolated bremsstrahlung amplitude 
and it leads to the characteristic GB radiation spectrum, H 2 + RBf . 

At finite ti and to — the numerical evaluation of the eq. (||) leads to the angular 
distribution of the radiated gluon. Fig. 2. displays our numerical results, normalizing the 
one gluon radiation one scattering distribution, dN^ /dyd 2 k±, by the one gluon radiation 

no scattering distribution, dN^ / dyd 2 k±. Fig. 2a, displays the contribution of the time- 
independent, classical cascading terms and their sum. Fig. 2b, shows the strong influence 
of the interference terms. Especially at small kj_ the hard radiation limit is recovered 
through the suppression of the GB contribution. This can be proved exactly in eq.([|) by 
means of k± — * 0. At large k± kinematical cuts appear from the appropriate jet energy, 
however the enhancement connecting the presence of the C\ cascading term can be seen 
clearly. 

— * 

In the case of more scattering (n s > 1) the existence of Cj cascading terms leads to 
approximate exponentiation of the relative gluon radiation distribution. Naive interpre- 
tation of this result would suggest that the energy loss exponentiates with the length. 
However, in that case multiple emission processes have to be also considered and this ne- 



cessitates a renormalization of the wavefunction as discussed in detail in Ref. 11 . With 



that renormalization one finds approximately constant dE/dx for n s < 3 up to L ~ 3 fm. 
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Fig. 2. (a) The normalized time-independent envelopes and separate contribu- 
tions to the conditional probability distribution of gluons associated with a single 
rescattering n s = 1 of quark jet (Eqj^j = lOOGeV) is shown as a function of the 
logarithmic "angle", log 10 A^_//i 2 with u = 0.5 GeV. These curves correspond to 
the incoherent parton cascade limit, (b) For finite t\o = 1 fm/c case destructive 
interference limits the corrections to the hard self-quench distribution to higher 
angles as shown for different gluon energies. 

Acknowledgment: We thank Urs Wiedemann for many discussions. This work was 
supported by the DOE Research Grant under Contract No. De-FG-02-92ER-40764, by 
the US-Hungarian Joint Fund No.652 and OTKA No.T025579, T024094. 

REFERENCES 

1. J. D. Bjorken, Fermilab-Pub-82/59-THY (1982) and erratum (unpublished). 

2. M. H. Thoma and M. Gyulassy, Nucl. Phys. B 351, 491 (1991). 

3. M. Gyulassy and M. Plumer, Phys. Lett. B 243, 432 (1990). 

4. M. Gyulassy, M. Plumer, M. H. Thoma, X.-N. Wang, Nucl. Phys. A 538, 37c (1992). 

5. X.-N. Wang and M. Gyulassy, Phys. Rev. Lett. 68, 1480 (1992). 

6. J. F. Gunion and G. Bertsch, Phys. Rev. D 25, 746 (1982). 

7. M. Gyulassy and X.-N. Wang, Nucl. Phys. B 420, 583 (1994). 

8. R. Baier, Yu.L. Dokshitzer, A.H. Mueller, S. Peigne, D. Scruff, Nucl. Phys. B483 
(1997) 291; ibid. B484 (1997) 265; ibid. B531 (1998) 403; Phys. Rev. C58 (1998) 
1706. 

9. B.G. Zhakharov, JETP Letters 63 (1996) 952; ibid. 65 (1997) 615. 

10. U.A. Wiedemann and M. Gyulassy, |hep-ph/9906257| . 

11. M. Gyulassy, P. Levai and I. Vitev, Jet Quenching in Thin Quark- Gluon Plasmas I: 
Formalism. 



